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1 Introduction 



Examining fundamental particle fields on the background of expanding universe, in particular de 
Sitter and anti de Sitter models, has a long history; special value of these geometries consists in 
their simplicity and high symmetry groups underlying them which makes us to believe in exis- 
tence of exact analytical treatment for some fundamental problems of classical and quantum field 
theory in these curved spaces. In particular, there exist special representations for fundamen- 
tal wave equations, Dirac's and Maxwell's, which are explicitly invariant under symmetry groups 
50(4.1) and 50(3.2) for these models. In the most of the literature, when dealing with a spin 
1 field in de Sitter models they use group theory approach. Many of the most important refer- 
ences are given below (it is not exhaustive bibliography, which should be enormous): Dirac [Il[2], 
Schrodinger (Sj H] , Lubanski-Rosenfeld [5] , Goto [6] , Ikeda [7] , Nachtmann [8] , Chernikov-Tagirov 
[9], Geheniau-SchomblondpiO] ■ Borner-Durr Tugov [T2], Fushchych-Krivsky |13] . Chevalier 
[E], Castagnino [Ml [19], Vidal [17], Adler [H], Schnirman-Oliveira pO], Tagirov [21], Rior- 
dan [22j, Pestov-Chernikov-Shavoxina [23], Candelas-Raine [24j, Schomblond-Spindel [23 l26] . 
Dowker-Critchley [27], Avis-Isham-Storey [21], Brugarino [22], Fang-Fronsdal [30], Angelopoulos 
et al [31] , Burges [32] , Deser-Nepomechie [33] , DuUemond-Beveren [34] , Gazeau [35] , Allen [36] , Fef- 
ferman-Graham [37], Flato-Fronsdal-Gazeau [38], Allen-Jacobson [39], Allen-Folacci [IQ], Sanchez 
[41j . Pathinayake-Vilenkin-Allen [42j, Gazeau-Hans ^43j, Bros-Gazeau-Moschella [44], Takook 
[15], Pol'shin [li iZl HI], Gazeau-Takook [l9], Takook [50], Deser-Waldron [SD [52], Spradlin- 
Strominger-Volovich [53], Cai-Myung-Zhang [54], Garidi-Huguet-Renaud [55], Rouhani-Takook 
[56] , Behroozi et al [57| , Huguet-Queva-Renaud [58] , Garidi et al [59] , Huguet-Queva-Renaud [60] , 
Dehghani et al |61j . Moradi- Rouhani-Takook [62], Faci et al |63] 

In a number of works, they have examined exact solutions of the wave equations on the back- 
ground of the de Sitter space and in particular the problem of passing quantum-mechanical par- 
ticles of various sorts (with different spin values and masses) through de Sitter horizon: Lohiya- 
Panchapakesan [69\ [70] , Khanal-Panchapakesan [71\ [72j , Khanal |73[ , Hawking-Page Otchik 
[75] . Motolla [76], Bogush-Otchik-Red'kov [77], Mishima-Nakayama [78], Polarski [79], Suzuki- 
Takasugi [80], Suzuki-Takasugi-Umetsu [STl [821183] 

In that studying, among other things, they had used a special definition for the so-called re- 
flection and transmission coefficients; then they had calculated corresponding quantities which in 
turn were used in considering the Hawking's formula [64l (Ml (Ml IST] as applied to the de Sitter 
space-time model 

exp(4|) - 1 ' ^ = ' ' 

where i?^ and F^- are reflection and transmission coefficients, respectively. Furthermore, on these 
lines, some rather complicated analytical expressions for those coefficients depending on the parti- 
cle's spin, mass, and quantum numbers (e, j) have been found. In contrast to most of considerations 
, the work by Otchik [75j has predicted (the case S = 1/2, m ^ was meant) that R^j = 0, F^ = 1. 
It should be noted that the outcome of the work [75], R^j = 0, has appeared accidental one in the 
sense that we have obtained the zero-value for a quantity that needs to be calculated. 

To date, any detailed analysis of possible explanation for that discrepancy has not been given 
yet. The present paper aims to investigate just those aspects of the problem. In the author's 
opinion, the problem to be solved is to show up the main points in conventional algorithm for 
finding those R^j and F^-, as applied to the de Sitter space-time, which provides subsequently 
the source of trouble. Such revealing study seems important because this general algorithm for 
calculating coefficients R and F was adopted (rather formally and without due caution) from 
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other and more complicated (and really barrier-contained) situations where any exact solutions 
are not known) : Staroninskiy - Churilov |841 [85] , Teukolsky - Press [86 , ST, 88j , Bardeen - Press 
[S2], Bardeen - Carter - Hawking [90], Unruh [OT], Fabbri [H^, Wald [QSj, Boulware [Hj, —age 
[951 |96l [97], Chandrasekhar - Detweiler [98], Matzner - Michael [99], Guven [10^, Bekenstein - 
Meisels [lOlj , Martellini - Treves |102j , Jyer - Kumar [103j , Hawking - Page |104j , Chandrasekhar 

Let us outline the content. In Sec. 2, we briefly reexamine early known solutions of scalar 
field in the de Sitter static coordinates' background. In so doing, certain their properties, the most 
material to the following, are itemized. In particular, four different solutions are singled out for 
special treatment (see Sec. 3) 

the notation used will be spelled out below. Here, / and ^f) are real-valued and linearly independent 
solutions of the matter radial equation (see (2.3)); respectively regular and singular ones at r = 0. 
The functions and again are linearly independent ones, complex- valued, and conjugated 
to each other. As a matter of fact, these pairs /(r) , g{r) and U^^{r) , U^.'^'^^r) can be related 
to each other by means of a linear transformation. So, there are some grounds for thinking that 
the above solutions ^ltand(^) ^ltand(^) represent standing waves, as well as ^>°^^ (x) and 
(x) describe running waves. Besides, as noted in Sec. 4, asymptotic behavior of all these 
solutions agrees with the used terminology. In addition, we will demonstrate that as the de Sitter 
curvature radius tends to infinity {p — )• oo), these de Sitter space-time model's solutions are reduced 
to the well known standing and running waves in the flat space-time background. 

In addition, as shown in Sec. 5, a radial component of the scalar particle's conserved current 
Jr{x) vanishes for the solutions ^^tand (^) ^ stand (^)' whereas to the (x) and (x) 
there correspond non-zero current waves; moreover, the relation {JrY"^' = —{Jr)°^^' holds. 

Finally, it should be added that a form itself of the effective potential curve, with no barrier, 
for the scalar particle's radial equation (this matter is treated in more detail in Sec. 6) might 
be able to refute absolutely the formulation itself about the need to study the process of passing 
quantum-mechanical particles through de Sitter horizon. 

So, the situation seems like as if the whole of the flat space-time model's case is carried, only 
with slight alterations, over de Sitter's model. Correspondingly, quite justified suspicions that some 
of the early established results on the process of particles' penetration through the de Sitter horizon 
are partly wrong might arise. However, it seems difficult to discern at once which part of the early 
performed study gives rise to subsequently misleading conclusions. In the same time, the mere 
claim that those results cannot be correct could be hardly enough. 

In sec. 7 we discuss how approximations in the form of wave functions can influence physical 
results. In sec. 8—10 we extend results to the case of electromagnetic field in de Sitter space. 
Finally, in sec. 11—12 we briefiy consider the problem for spin 1/2 field in de Sitter space. 

The main goal of the present paper is to clarify the whole situation with the particle penetration 
through de Sitter horizon. The paper is partly a pedagogical review, so we intend to describe in 
detail the mathematics used to treat particles in presence of curved space-time background. Our 
treatment contains much technical details that may be of small interest for skillful readers but 
helpful for beginners. 
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2 Solutions of a radial equation 



Conformally invariant wave equation for a massive scalar field ^{x) in de Sitter space has a form 
(M = mcp/h, p is the curvature radius): 



We shall analyze solutions of that equation in the static coordinates 

dS"^ = ^dt^ ^^(^"^^ + ^^'?^^) > < ^ < 1 > ^ 
Taking ^{x) in the form of a spherical wave 

^{x) = e-^^* /(r) Yj^{e, </.) , e = Ep/hc 

for /(r) we get 



(2) 



(3) 



d^f ,2 ^>',d/ ( M2 + 2 + 



$2 



/ = o. 



(4) 



All solutions of that equation dH) can be expressed in terms of hypergeometric functions. To this 
end, after introducing a new variable r"^ = z, the equation reads as 



d^ d 
Hl-z)z^ + (6-10z)f + 

oz^ az 1 — z 



j{j + 1) 



/ = o. 



(5) 



and further taking the substitution f{z) = z'^ {1 — zY F{z) , where parameters {k, a) are to be 

K = j/2, -(i + l)/2, a = ±ie/2, 

for S{z) we get 

z{l - z)F" + [c-{a + h + l)z\ F' - ab F = 

where (a, b, c) satisfy 

c = 2k + 3/2, a + b = 2K + 2a + 3/2 , 

2 2 3 3 M2 + 4 

ab = K + 2k(T + (T H — K -\ — a -\ . 

2 2 2 

Taking k = j/2 and a = — ie/2, for (a, 6, c) we get a regular (at r = 0) solution 
AC = j/2, a = -ie/2 ,c = j + 3/2, f{z) = z^/^ (1 - z)-''/^ F{a, b, c; z) , 



3/2 + j + i\/M2-l/4-ie , 3/2 + j - iv'M2 - 1/4 - ie 

a = , = . 

2 ' 2 

In turn, choosing k = —(j + l)/2 and a = —ie/2), we will obtain a singular (at r = 0) solution 

K = -{j + l)/2 , = -ie/2 , c=-j + 1/2 , 
1/2- j + iy'M'^ -1/ A- ie 



(6) 



a 



/3 



1/2 - j - «VM2 - 1/4 -ie 



;(^) = ^-a + l)/2 (1 _ ^)-«/2 ^(^^ ^. ^) _ 



(7) 
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3 Standing and running waves 

It is easily verified that the regular solution above, first, represents a real-valued function; second, 
it can be resolved into two complex-valued (and conjugate to one another) solutions of the same 
equation (jj]). Applying here the terminology used in the ordinary case of spherical waves in the 
flat space-time model, one might say that the regular standing wave is certain superposition of two 
running waves; one of the latter travels to the de Sitter horizon (out-wave), another runs backwards 
it (in- wave). The same is valid for the singular standing wave: the g{r) is a real- valued alike, and 
it can be resolved into a linear combination of the same complex-valued solutions (certainly with 
other coefficients). 

Now, let us examine this matter in full detail. First we shall consider the regular solution. To 
realize a required decomposition, it is sufficient to use one of the so-called Kummer's relationships 



r(c)r(c-Q-b) r(c)r(Q + fe-c) 



T{c-a)r{c-b) 



T{a)T{b) 



(8) 



where Ui, U2, Uq represents three different solution of the same hypergeometric equation; those 
three are determined by 



Ui = F{a, b, c; z) , U2 = F{a, b,a + b - c + l;l - z) , 
Uq = F{c — a,c — b,c — a — b + 1;1 — z) . 



(9) 



Applying that relation 

/(^) = 

where 



dH]) to the solution f{z) as Ui, we shall get 

r(c)r(c - a - 6) , . . r(c)r(a + 6-0 



T{c-a)T{c-b) ' 



r(a)r(6) 



jjtn. 
^ run. 



(10) 



TjOUt. / N 



z-"- (1 



J7*" (z) 



J/2 
z3l^ (1-Z 



- z)-"/2 F(^a, b,a + b-c + l;l- z) , 
F(^c - a,c-b,c - a-b + l;l- z) 



(11) 



Noting that 



(a - a) , b* = (c-b) , {a + b- c)* 



(a + b 



one can conclude that those two functions U°^^{z) and U^un.i^) referred to each other by 
the operation of complex conjugation 



(^rwn.(^)) ~ ^rwn.(^) ■ 

Besides,the function f{z) being a sum of two complex-conjugate expressions (see (fTO|) and p^ ) is 
a real-valued one. Hence, one can get 



fiz) = 2 Re 



r(c)r(c -a-b) 
r(c-a)r(c-6) 



Uout{z) 



2Re 



T{c)T{a + b-c) 

mm 



Uiniz) 



(12) 



To obtain an analogous expansion for the singular function g{z) (2.6) , one ought to apply the 
same Kummer's formula dH]), but now having replaced (a, 6, c) by (a, (3, 7): 



r(7)r(7-a-/3) r(7)r(a + /3-7) 

Ul = —, T^TT -r- U2 + Ufi 



r(7 - a)r(7 - /3) 



r(a)r(/3) 



(13) 
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where 



Ui = F{a, (3, 7; z) , U2 = F(a, /?, a + /3 - 7 + 1; 1 - z) , 

i76 = F(7-a,7-/3,7-a-/3 + l;l-z) . (14) 



Thus we arrive at 



where 



_ r(7)r(7 -a-/3) , r(7)r(a + /3 - 7) 

^^^^ " r(7 - a)r(7 - /3) ™" ^ ^ ^ r(a)r(/3) ^'■'^'^-^^^ ' ^^^^ 



f^™n.(^) = z'^' (1 - - a, 1 - /3,7 + 1 - a - /3; 1 - z) (16) 

and the following relations 

a + 1 — 7 = a, /3 + 1 — 7 = 6, a + ^ — 7 + l = a + 6 + l — c, 

(1-a) = (c-a) , (1-/3) = (c-6) 

are taken into account. Again, the function g{z) is real-valued, and one can get 



g{z) = 2Re 



r(7)r(7 -a- /3) ^ 



r(7 - a)r(7 - /3) 



2Re 



r(7)r(a + /3-7) , . 

r(a)r(/3) ^run.yz) 



(17) 



4 Asymptotic behavior and fiat space-time limit 



Now let us consider how the above four solutions behave at two natural asymptotics: at the origin 
r ~ and at the horizon's region r ~ 1. First, let us single out for that treatment: 

UZniz) = z^'"^ (1 - z)-''/'^ F{a, b,a + b-c+l;l-z) . (18) 

Again, applying the formula ([8]) to U°^^ {z) as Ui, we arrive at the following expansion 

(z) = r(a + fc + i-c)r(i-c) /2 _ ,,/2 

t^rnnA^J r(5 _ c + i)r(a - c + 1) ^ ^ ^{a,o,c,z) 
+ "^^"^'pA^rr'"'^ ^^"■"'^'^ (1 - ^)"^^'' F(6 - c + 1, a - c + 1, -c + 2; .) . (19) 

r(a)r(6) 

that is appropriate for asymptotic analysis: 

U^t (r ~ 0) ~ , UZL (r ~ 1) ~ (1 - r2)-/2 . (20) 
Introducing a new radial variable r*: 

* Pi 1+^' exp(2rVp)-l 

r = - In , r = -. , r e 0, oo) 

2 1 - r ' exp(2r*/p) + 1 ' ^ ' ^ 

one can easily reexpress the second asymptotic relation in (j20p in the form (e = Ep/fic) 

UZnir* ~ oo) ~ f2-iEp/nc\ ^^p^+iEr* /he) . (21) 
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Analogously, one can find asymptotics for the in-wave: 



f3 



or 



UZn.ir* ~ oo) ~ (2+^Ep/nc\ ^^p^^_iEr* /he) . 



Proceeding with that analysis for the standing waves, we can easily get 

1 



f{r ~ 1) 



f{r ~ 0) ~ r-* 

r(c)r(c -a-b) 
r(c-a)r(c-6) 



g{r ~ 0) 
(1 - r^)— /2 



+ 



r(c)r(a + 6 - c) 



(22) 



(23) 



(24) 



2 i?e 



r(a)r(6) 

r(c)r(a + 6 - c) 



:i -r 



exp{—iEr* /he) 



g{r ~ 1) 



r(a)r(6) 

r(7)r(7 - « - /3) . 2wi./2 , r(7)r(a + /3-7) _ W2 

^ ^ rfa)r(/3) ^ ' 



(25) 



r(7 - a)r(7 - /?) 

= 2 i?e 



X«)r(/3) 

r(7)r(a + /3-7) , .^^^ ^ 

r(a)r(/3) ^ """"P^ '^"^ /''''^ 



(26) 



Now, we take up examination of a flat space-time limit as applied to the de Sitter's solutions. 
For definiteness we single out the above out-wave. It will convenient to start with its representation 
according to (|19|) . At this, one essential step should have been done in advance; namely, one ought 
to have rewritten the above-used expressions for parameters (a, b, e) with a curvature radius p 
apparently separated and set aside (a subsequent limiting procedure will be realized just through 
that parameter: p ^ oo ) 



p+l-iep + ivV^Vf2^^l74 



p+l-iep- iy/p^M^ - 1/4 



(a - e+1) 



e+1) 



-p+l-iep + iy/p^M^ - 1/4 



-p+l-iep- i^^p'^M'^ - 1/4 



p = j + l/2 ; (a + fe + c) = (l-iep) , 

and limp_j.oo(/0^-z) = , where R is the usual radial coordinate of the flat space-time. 
Using a conventional series representation for a hyper geometric function 



(27) 



F{a, b, c; z) 



ab z a(a + l)b(b + 1) z^ 

1 + -- + ^ ^ + 



c 1! ' c(c + 1) 
and further taking into account two limiting relations 



(UI<1) 



(28) 



lim 



(a + n)(fe + n) ^ 2 
(a-c + l)(b-c+l) 2 2 

_9 P 



lim 



_(e2_^2^ _ 



1 



.(e'^-M^) - 



7 



one comes to ((e — M^) = k^): 



\im F{a,b,c;z) =T{l+p) Y] 
p->-oo ^n\L[l + n + p) 

(-A;VV4)" 





oo 



Urn F{b- c+l,a- c + 1, -c + 2; z) = r(l - p) V 



p— i>00 



n!r(l + n — p) 



(29) 
(30) 



Now, with the use of a known expansion for the Bessel functions 



from (I29])-(I30D one gets 



hm 

p—>-oo 



pi (1 _ ^2)-^./2 ^(^^ ^. ^) 1 ^ ( ^ ^) 

J k 

J-p{kr) 



hm 



p 



(31) 



Further, taking in mind ([27]) and ([31]) -([HI]), for U°^^ , at the hmit p — 7- oo, we obtain {A is yet a 
non-fixed constant) 

hm A U^^ {z) = hm ^ X 

p— >oo ■ p— >oo 

r(-k/) + 1) r(-p) r(i + p) (2/fc)p jpjkr) p-p+y^ i_ 

rlli+i^/p'^M^ - 1/4 -iep + p + 1)] r[i(- Vp^^^ -l/i-iep + p+l)] 7^ 
r{-iep + 1) r{+p) r(l - p) (2/A;)-P J-pikr) p+v+^l^ 1 

r[i(+iv^p2^2 _ _ _ ^ + 1)] r[i(-i^p2^/2 _ _ _ p + i)j ^ 

Now, using the known relation for F-function 

T{p) F(l-p) = -^, 

and introducing shortening notation 



-iep + iJp^M"^ - 1/4 M-e 

1 ~ 2 



-iep-iJp^W^ - 1/4 M + e 
P^i ) 7. ~ ^ — 7, — 



PZ2 



(32) 



we change (j32]l into 



vr T{—iep + 1) 



hm yl 
p-i>oo sin(7rp) ^Jr 



hm A (^) 

p— ^-oo 



(2/A;)P Jp{kr) p-P+^l'^ 



— J—— J— + (P ^ -P) 

F(/92;i - T{pz2 - 



(33) 



the designation [p — )• —p) means that here else one term is to be placed, which is derived from the 
previous one with a change p — t- —p. The constant A will be chosen to guarantee the existence of 
a finite limit in (4.9); so it is 



p-j-oo p-s-oo F(a + o — c+1) p^-oo F(— iep + 1) 



(34) 
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One may carry out further calculations solely with the first term in (|33p (it will be designated as 
F{p) ; the second term's contribution will be determined by the first through a change p — t- —p) 



F[p) = lim 

p— ^oo 



-Ti 1 {2/k)P Jp{kr) r{pzi + 1/4) r(pz2 + 1/4) 



sin(7rp) ^ 2P-1/2 Y{pzi - ^) T{pz2 - ^) 



(35) 



So, let us dwell on the term F{p). Applying else one time the above formula referring F-s with 
arguments p and (1 — p), one can produce 



F(p) = lim 

p— >-oo 



-vr 1 {2/k)P Jp{kr) r(/9Zi + l/4) sin7r(pz2 + ^) T{-pz2 + ^) 



sin(7rp) ^ 2P-1/2 r(pzi + i^) sin7r(pz2 + 1/4) T{-pz2 + 3/A) 
which after having used the known asymptotic relationship for F-functions 



(36) 



oo 



argx < TT 



F(a; + a) 
F(x + /3) 



and taken into account the equality ((e + M) > 0) 



smTT{pz2 + 



1— p ^ 
2 ^ 



exp [vr 



p(M+e) 



+ ivr 



i+Pi 



exp[— TT ^^^^"'"'^-* 



~ X 



ITT 



sin7r(/9z + 1/4) 
provides us with the expansion 



exp[7r^^^^ + m/4] - exp[-7r^^^^ 



m/A\ 



p i 

lim exp[i7r(- + -] 



F(p) 



lim 

p— >oo 



-vr 



1 



^ ^ {^Y Jp{kr) exp[(pzi)^'/2-V4 (p^2)^/'-^/1 exp[ivr(f + \] 



from where it follows 



V fcr 2 2 smvrp 



(37) 



Remembering the second term's contribution {p — )• —p), the final limiting expression for A Uoutiz) 
will be 



lim A Uoutiz) 

p— >oo 



(38) 



where Hj^^-^^^^{kr) denotes the Hankel function 



^j+l/2(^) 



ip 



sin(7rp) 



[ e'P^ Jp{x) - J_p(x) ] 



(39) 



The right-hand part of the relation (j38p provides quite standard representation for an expanding 
spherical wave (going to the spatial infinity) in fiat space-time model. Thus, the notation in the 
Sitter space model 



^sing. (h'^f'g- (f,out. 

stand. ' stand. ' ^run. ; ^run. 



is correct and consistent with that used in the fiat Minkowski space. 

In connection with that limiting passage from the curved to fiat space-time, one point of no- 
ticeable importance for the following should be noted. Because of the curved metric ^ practi- 
cally coincides with the flat (Minkowski's) one as soon as one gets far enough from the horizon 
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(0 < r « p), one might expect little difference between '"curved" and "flat" solutions of 
the respective (scalar particle's) equations: fcurve{i") ^ ffiat{f) ■ Being applied to a solution of the 
kind of extending ones, that claim can be formulated as follows: the approximation 



1 y kr -J+i/s*^ ' 



(40) 



is quite valid. However, as will be seen, such a suggestion is not completely right: some limitations 
accompanying (j40p will be readily discerned. 

To clear up this matter, we shall compare carefully two respective equations: 



in the de Sitter model , 

d^ 2(1 -2r2) d e2 

+ — ^ — + 

dr^ r(l — r^) dr (1 — r^)^ 

in the Minkowski model , 

2d_ . 
dr"^ r dr 



^ = e-''' f,,{r)Yj^{e,. 



M2 + 2 j{j + 1) 



1 — 

' - M 



itt M 
2 j{j + 1 



(41) 



At the region r << p , the equation (|4ip may be rewritten as follows 



d^ ^ 2 d 
dr'^ r dr 



r , o , 2r^ 

- + ...+ 1 + — + .. . 



(m2 + -)(1 + - + ...) 



j{j + 1) 



(1 + - 



■) 



. 



(42) 



On comparing (I42p with ()4ip one may conclude that those practically coincide if and only if two 
inequalities 



r « p and | 
hold; from the latter it follows 



(43) 



(44) 



That is, the difference between the particle's energy e and its mass m should be great enough as 
compared with a scale given by (j/p), where p denotes the curvature radius, and j is an orbital 
momentum number. At j = 0, the relation (e^ - M^) >> 2/p'^ ought to be held. The curvature 
parameter p is quite a big but finite one, whereas no limitations on j exist. 

5 Standing, running waves and conserved current Ja{x) 



As an additional argument to justify the above-used terminology let us consider properties of 
the scalar particle's conserved current 



J„(x) = i(^>*V„$ - ^*Va^) 



(45) 
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Using a general spherical wave's representation <I>(x) = e f{r)Yjfn{d,4') ) foi' current's compo- 
nents Ja{x) one gets 

Ji = 2e I / |2 I |2 , J,^ = -2m | / p | Gj™ \^ , 



Jr = iir-^f - f^n I 0jm P , Je = (46) 



whence it follows that standing and running solutions (in the above terminology) are characterized, 
respectively, by 

( J Y'^a- = Q (J ying. ^ ^ , J .out. - _ (J \in. 

stand. — ^' y'^^'lstand. — ' K'^rJrun. — \'Jr)run. ■ 

Those properties of a radial current component are in agreement with our notation which prescribes 
consider $r^J;f*""-*-solutions as relating to radial in- and oui-streams of particles, whereas ^Itand^^^ 
represent radially immobile particles' states. It seems likely that analogous relationships conforming 
to the anticipated behavior of such waves will be established if we turn to the energy-momentum 
tensor referring to those solutions. Thus, the terminology exploited above seems quite reasonable 
and appropriate. 



6 On transmittivity and reflectivity of the de Sitter horizon 



From the preceding material it might be concluded that the phenomenon itself of any reflection 
of scalar particles by the Sitter event horizon should not be observed whatsoever. Added credence 
to that opinion may be given by qualitative but quite heuristic consideration of the problem in 
terms of an one-dimensional Schrodinger's like equation with an effective potential U{r). Indeed, 
the matter radial equation @ can be taken into the form (the variable r* is used again) 



dr 



*2 



U{r*)) 



G{r*) = ; 



(47) 



where G{r*) is given by 

/(r) = ^l{r*) G{r*) , 



dr* 



In ^ 



(1 -r)(l + 2r) 



P 



and an effective potential U{r*) is determined according to 



U{r*) 



1 



4(1 -r) + 



1 + r 



2 2 , J(i + 1) 
+ m p + K 



(48) 



As readily verified, that potential (j48p creates an effective repulsive force —dU/dr* at every spatial 
point (r < 1; r* G [0, +oo]) of the de Sitter space-time. On other words, the center r = effectively 
repulses the scalar particle everywhere and no barrier appears between the center and event horizon 
region r ~ /j: 



dU 
dr* 



1 



+ 



1 + r 



+ 4(1 -r) + (1 



+ 1) 



2r 
+ 4 



j{j + 1) , 2 2 

5 + m p 



1 



(1+r) 



> . 
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In a region near to the horizon (r* ~ +00), the potential U{r*) tends to zero, and correspondingly 
solutions (running waves are meant here for definiteness) look like exponents 

G{r*) ~ exp(ibzer*) ; 

the latter boundary behavior may be naturally interpreted as associated with waves passing to (+) 
and from (— ) the de Sitter horizon (r ~ 1 or r* ~ 00). 

After having examined this Scrodinger's-like equation, it has become absolutely evident that the 
quantity such as a reflection coefficient R^j cannot be correctly determined. However, it would not 
be sufficiently enough if only this single statement had been formulated. In so doing, undoubtedly, a 
bit of vagueness would still remain. Because of this, now let us reexamine closely (if not narrowly) 
an algorithm for calculating such a quantity R^j on the background of the de Sitter space-time 
model. For definiteness, let us consider the case of massless particles. 

A wave going to the horizon is described by the following radial function (see (fT9ll at = O 
and r = R/ p) 



T(a + b+l-c)T(l 



T{b- c+l)T{a- c+1 



+ 



r(a + 6+ 1 -c) r(c- 1) 1 



r(a) r(6) 



Y (1 - r^)-"/2 F{b- c+l,a- c+l,-c + 2;r'' 



(49) 



where 



J - 



i + 1 - iep 



j + 3/2=p+l. 



At great distance away from the horizon (R « p), for the (i?) one has an approximated 
representation 



^run. (^) 



r(a + 6+l-c)r(l-c) 1 2 . JpjeR) 

T{b-c+l)T{a-c+l) pi^e' ^^^^P) 



R 



r{a + b + l-c) r(c-l) ,+1 2 _ J_p{eR) 

r(a)r(6) ^ ^e' ^ 



(50) 



but one ought to remember about accompanying limitation (j44p on quantum numbers (see above 
the limiting passage from the curved to flat space-time model): 



R«p, j «ep. 
This means that to be correct one will must take the limits 



lim 

ep»j 



T{a + b + l-c)T{l-c) 
T{b-c+l)T{a -c + 1) 



lim 

ep»j 



r{a + b + l- c)r(c- 1) 

mm 



(51) 



(52) 



if one is going to exploit that expansion (I50p subsequently in a noticeably significant physical 
context. For the moment, setting away the need to realize additionally the above limiting passages, 
let us proceed with calculations providing "a required" analytical expression for R^j, following 
the standard pattern. So, taking into consideration an asymptotic formula for Bessel functions 



X » z/" 



Jix) 



r(2;^ + l) 2-^"-^/^ J_ 
T{u+1) r{iy+l/2) 7^ 



exp ( +i{x - ^{u + ^)) ) + exp 



-i[x 



(53) 
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we get 



where A and B are 



\-ieR 



eR 

leR 



j < J « eR « ep 
1. 



A exp(-z|(p+ ^)) + B exp(-«| 



-p + 



eR 



A exp(+i|(p+^)) + B exp(+i^(-p + ^)) 



(54) 



B 



_ r(a + 6 + 1 - c) r(l - c) 2-^-^ T{2p + 1) 
~ r(6-c+l) r(a-c+l) (e/9)i r(p+l/2) 
r(a + 6 + 1 - c) r(c - 1) (epp+i r(-2p + 1) 



r(a) r(6) 



2-j rip + 1/2) 



(55) 



The reflection coefficient R^j, by definition, is a square modulus of the ratio of an amplitude at 



e "^^/ei? to an amplitude at e'^ /eR ; thus we shall find 



R. 



I e 



-nrp 



[A/B) e^ 



■?.7rp 



+ 1 



[A/B) e-'-^P + 1 ■ 

It remains to show that, after having into account the limitation ep >> j, this expression (|56p for 
will identically vanish. For one has 

A _ 1 r(a) r(b) r(i - c) r(-p + i/2) r(2p + 1) 

~ (2ep)2i+i r(a-c+l) r(6-c+l) r(c - 1) r(p+ 1/2) r(-2p + 1) 
After using the formulas 

T{p) T{l-p) 

we find 



vr 



r(2x) 

Tlx) 



sin(7rp) 

r(i-c) r(-p + i/2) r(2p + i; 



r(z + i/2) 



r(i/2) 



(57) 



(58) 



r(c-i) r(p + i/2) r( 

Further, applying the asymptotic relation 



we get 



lim 

ep»j 



{ X — )• DO , I arg X I < TT } 

r(a) r(6) 



-2p + 1) 
r(x + a) 

^ r(x + /3) 



_24i+2 _ 



(59) 



^ 2 ' 



\2p 



r(a-c + l) r(6-c+l) 
Substituting ([59]) and ([60]) into (j57|) . for ^/i? we finally arrive at a simple result 



lim 

ep»j B 



-i-i) 



2p 



(60) 



(61) 



and therefore the expected identical zero for R^j arises, R^j = 0. 

In other words, analytical expressions for R^j existing in the literature are due to not taking 
the all required limits in the final expressions. It should be stressed that, as can be seen from the 
present study, certain analogous limiting procedures as well are to be performed in other and more 
complicated physical cases. Correspondingly, the possible underestimating of all set of limiting 
conditions in every separate situation , coupled with calculation of final analytical expressions of 
coefficients R and P, might be source for partly incorrect conclusions. 
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7 Approximation influencing physical results 



To prevent anybody from possible errors let us discuss the formula ()50p giving decomposition for 
a wave going to de Sitter horizon at far distant domain and consider possible discrepancies arising 
from the use of incorrect approximate expressions for wave functions. Indeed, instead of correct 
decomposition (|50p let us write its (by hand) modification through introducing at Jp{eR) and 
J-^p{eR) some coefficients (1 + Ap) and (1 + A_,p) (to correct formula there correspond Ap = and 
A_p = 0): 



C/o«t(i?)~ [(1 + A, 
+ (1 + A_p) 



T{a + b+l - c) r(l - c) 12 



r(6-c+l)r(a-c+l) pi \ 

2 
e 



. (-rr(i+p) -^^^'^^ 



R 



T{a + b+l-c) r(c-l) 

mm ^ 



r(i-p) 



J-p{eR) 



R 



Correspondingly, instead of A and B according to (I55p we have 

bA ^ A' = {1 + Ap) A , B ^ B' = {1 + A^ 



B 



and further 



nm — —7- nm — 

<^p»3 B <^p»j B 



lim 



1 + A, 



ep»j 1 + A_p 

Therefore, now we arrive at a new expression for reflection coefficient (compare with ([56 

{A'/B')i+i)^P + l ,2 



(62) 



(63) 



(64) 



R',j 



{A'/B'){-i)2p + 1 



Ap - A_p 
2 + Ap + A_p 



(65) 



In other words, just an error in approximate form for a wave function leads us to a quite "physical 
result" Such incorrectness may be presented if any uses approximate solutions for wave eauations 



R ~ lim 



ifP + 1 



(66) 



8 On Maxwell equations in general covariant tetrad Majorana — 
Oppenheimer form on the background of de Sitter space 

Below we show that the main results can be extended to the case of electromagnetic field. At 
this we will use an old and almost unusable in the literature approach by Riemann -~ Silberstein - 
Majorana - Oppenheimer in general covariant tetrad form. Maxwell equations in Riemann space 
can be presented as one matrix equation |1061 11071 HOSj 



a 



-il 





E + icB 



^ = J{x) 



a 



10 
-10 

1 



, a 





-1 




1 

1 



-10 



, a 



J = 


1 


P 




eo 


ij 








1 





-1 





1 








1 









(67) 
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or 



^(e?o)^P + l^f'labo)^ + aHe'^,)dp + -f'jabk)^ = J{x) . 



Allowing for identities 
1 



[si(7230 + «7oio) + •52(7310 + «7o2o) + 53(7120 + no3o)] 



1 



f^lahk = [Sl(723fc + «701fe) + S2(731fe + «702fe) + S3(7l2fc + ^703fc)] 



where 









Sl = 




n 


,S2 










n 


, S3 = 







n 





















1 









-1 





n = 








-1 


, T2 = 











, T-3 = 


1 













1 







-1 






















and using the notation 



(0)^P - 9(0) , 



do = <9i 



(fc) 



a = 0,1,2,3 , 



(723a, 731a, 7l2a) = Pa 



(701a,702a,703a) = Vq 

eq. (f68l) in absence of sources reduces to 

- i[ 9(0) + s(po + ivo) + [ 9(fc) + s(pfc + iv^) ] ^ = 
Let us consider this equation in the de Sitter static metric and tetrad 

dS^ = $ (it^ - — - r2(d6'2 + sin^ Od^'^) , $ = 1 - , 

efo) = (^,0,0,0) , ef3) = (0,V¥,0,0), 



^(1) 



(0,0, -,0) 
r 



'(2) 



(1,0,0, 



1 



7030 



7311 



2\/$ r 
we get to explicit form of a matrix equation 



7322 



r sm t) 

^ cos 6 

, 7122 = — ^ 
r smt^ 



+ ^f^{o?dr + + —S3 + -Se,^ 







a 



^e.d = — de + a" 



ds + S3 cos 6 



sm ( 



(68) 



(69) 



(70) 
(71) 



(72) 



(73) 



It is convenient to have the spin matrix S3 as diagonal, which is reached by simple linear transfor- 
mation to the known cyclic basis 



U 



-1/^/2 i/V2 
1 

l/^/2 



Ua 



1 

u 



-1/^/2 l/^/2 
-i/\/2 
10 



(74) 
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so that 



1 

71 





1 

V2 





— z 







1 

"71 





-1 











+1 








r[ = 


— i 





—i 


—i 





1 






1 


-1 






— i 














-1 





-1 





1 









—i 





—i 










1 





1 





—i 





'2 


1 


—i 





-1 





,a ^ = 





—i 











— i 





— i 


"71 





1 





-1 


-1 











-1 





— z 







— i 





1 
















+i 



Eq. ([73]) becomes 



'2 J 



+ \/$(a ^dr + ^ + 



2^> 







a 
r 



'1 



+ a ■ 



'9 5^ + s'-i COS ( 



S1I1( 



(75) 



9 Separating the variables and Wigner functions 

Spherical waves with (j, m) quantum numbers should be constructed as follows 



-iujt 





f2ir)Do 
/3(r)D+i 



(76) 



where the shorted notation for Wigner D-functions |112] I113j is used: D„ = Di 



0) , a 



— 1,0, +1; j,m determine total angular moment. We adhere notation developed in [12^ I12I1 \1'22\ 
I123j ; before similar techniques was applied by Dray I117j , Krolikowski and Turski [118j , Turski 
[119j ; many years ago such a tetrad basis was used by Schrodinger [ 109j and Pauli [llOj when 
looking at the problem of single- valuedness of wave functions in quantum theory. In the literature 
equivalent techniques of spin- weighted harmonics Goldberg et al |114j (se also in |115j ) preferably 
is used though equivalence of both approach is known [116] I117j . 

With the use the recursive relations |113j (below u = j{j + 1) , a = sj {j — + 2)) 



m — cos ( 



sin( 



deDo 



1 Tft 

2 smu 

1 , , m + cos 9 
deD+i = -{vDq - aD+2) , ^ — 

2 sm& 



= -{aD^2 + T^Do) 
■- i {uDo + aD+2) 



(77) 



we get (the factor e is omitted) 



V 

72 



(/i + h)D^ 
-i f2D-i 
i (/i - h)Do 

+t f2D+i 



(78) 
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Turning to Maxwell equation (j75|) . we arrive at the radial system 

(1) V¥(-^ + -)/2 + -^(/l+/3) = 0, 

ar r r ^2 

(3) -^/2 + -^(/l-/3)=0, 

V<I' ^ v2 



(4) ( " +.v¥-;^ + ,:^ + i * )/3 + i-^/. = 0, (79) 



Combining equations (2) and (4), instead of ()79p we get 



(2) + (4) , -^(h + A) - + ^ + ^)(/, - A) = , 



-;iA+^^(/,-/3) = o, 

(1) + 2)A + + A) = , 

ar r r y 2 

It is easily verified that equation (1) is an identity when allowing for remaining ones. So independent 
equations are 

-^/2 + -^(/l-/3)=0, 

-^(A + A)-i(v/5;;^ + ^ + -2^)(A-A) = 0, 

$ dr r 2v^ 



-^(/, - A) - . (v/¥;f + ^ + ^)(/, + A) _ . ^ (80) 



Let us introduce new functions: 



then eqs. pOj) read 



p _ /i + /3 _ /i ~ /a 

^" ^/2 ' ^" V2 



/2 = 5, -^/-H:r + ~ + ^)f = 

w r <J> ar r 2<P 



The system (jSTj) is simplified by substitutions 



g{r) = G{r) , f{r) = F{r) 
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and it gives 



/2 = -^G(r), iu;F = ^^G, 



dr 



G = , 

So we have arrived at a single differential equation for G{r): 

£G dG fuj^ + l 



dr^ ~^ ^ dr ~^ \ $2 



G = 0. 



(82) 



(83) 



or taking $ = 1 — 



d^G 



2r dG 



+ 



jU + 1) 



G = 



(84) 



This equation coincides with that following from dH), if one translates eq. (j3]) at M = to new 
function /(r) = r^^Lp{r): 



d^ 



2r dip 



OJ 



3{j + 1) 



(1^.2 1 — dr \(1 — r2)2 ^2^]^ _ j.2^ 
With the variable z = r^, eq. ([85]) gives 

i(j + 1 



ip = Q. 



4.(1-.)— +2(1 -3.)- + 



G = 



(85) 



(86) 



and after substitution G = z"'{l — z)^F{z) we arrive at 



+ 



4.(1-.)^ +4 

4o2 -2a- j{j + 1) 462 _^ 



2a + ^-(2a + 26+^)z 



+ 



a; 



l-z 



4(a + 6)(a + 6 + 



1. 



dz 



F = 0. 



Requiring 



a = + 



i + 1 



2 ' 



2 ' 



we get 



2a + ^ - (2a + 26+ ^)z 



— -(a + 6)(a + 6+-)F = 0, 



parameters of hypergeometric function are given by 



1 1 

a = a + b , i3 = a + b + - , j = 2a + - 



(87) 



(88) 



(89) 



Evidently, these solutions coincide with those described in Sec. 2. 
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10 Spin 1/2 particle 



Dirac equation (the notation according |108j is used) 

i7'(efe)5a + ^CT^Sabc) - M] * = 

in static coordinates and tetrad of the the Sitter space takes the form 



^0 _ 1 31 , 2 -32 ^/ 



1 



S^^-M]^(x)=0, 



id + i a 



12 



cos ( 



sin( 



Below the spinor basis will be used 

/ 
/ 



-CTj 

a, 



1(7 



-12 



fT3 
^3 



Allowing for ^^a"^^ + 7^i^^ = 7^, 7°(T°^ = 7^/2, eq. ([9T]) reads 

^^7i^* +.^/l7=^(a. + i + ^) + is^,<^-M}M/(x) = 0. 
One can simplify the task with the help of substitution ^(x) = r^^<l>^^/^ F{x): 



( i + i\/¥73 a, + - - M ) F(x) = 



Spherical waves are constructed trough substitution D-L^ai^'i 0) = D^): 



-let 



Mr) D_y^ 
f2{r) -D+i/2 
fsir) D_y2 
/4(r) -D+i/2 



With the use of recursive relations |113] ) 



de 



1/2 



1/2 



a D_i/2 - b -D+3/2 , 



b D_s/2 - a -D+i/2 



-m — 1/2 cos 6 

sin0 
-m + 1/2 cos 9 

sin 9 



D 



+1/2 



D_ 



1/2 



-a i:'_i/2 - b 15+3/2 , 



-6 L'_3/2 - a -D+i/2 



where a = (j + l)/2 , b = (1/2) ^(j " l/2)(j + 3/2) , we get 



- Uir) D__y2 
+ fsir) D+y2 
+ /2(r) D_,/2 

- Mr) 15+1/2 



u = {j + 1/2) 
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And further we arrive at radial system 
d 



ar r ar r 



-^/i + iV^:rf^ + '-f^ - M/3 = , - iV^-^f2 - i-h - M/4 = o . 

ar r y$ ar r 



(96) 



To simplify the system let us diagonalize P-operator. In Cartesian basis it is lie. = ^S" P, after 
translating to spherical tetrad it looks 



n 



sph. 












-1 








-1 








-1 








1 

1— I 












p . 



(97) 



-iet 



From the equation Ilsph. ^jm = H ^ jm it follows H = 5 (— 1)-'"''^, 5 = ±1 and 

/4 = (5 /l, h = S f2 , '^{x)ejni5 = 

Allowing for (j98|) . we simplify the system (|96p : 



fi{r) Z)_i/2 
/2(r) 15+1/2 
'5 /2(r) D-i/2 
^fiir) D+i/2 



(98) 



+ - )/ + (^ + 5M)5 = 0, 
ar r v$ 

4 = 



(99) 



where instead of fi and /2 new functions are used 



/ 



/1 + /2 
x/2 ' 



5 



/1-/2 
i\/2 



For definiteness, let us consider eqs. ([99]) at = +1 (formally to the second case 6 = —1 there 
corresponds the change M =^ —M): 



{^j- + - )/ + (^ + M)9 = 0, 
ar r 



- - ) 5 - (^ - M) / = . 
ar r v/$ 



(100) 



Here there arise additional singularities at the points e + M = , e — \/5 M = 0. Correspond- 
ingly, equation for /(r) has the form 



dr"^ 



2r 



Mr 



l-r2 ^/rTr72(g + ^/^/i^); dr' V(i 



,,2^2 



+ 



Mu 



v{v + l) V 
'r2(l-r2) ~ (l-r2) ^1^72 ' ^1^72 (e + M\/r^^) 



M2 
1 — r2 

/ = o. 



(101) 
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However, there exists possibility to move these singularities away through special transformation 
upon the functions f{r),g{r) (see [75]). To this end, as a first step, let us introduce a new variable 
r = sin p, eqs. pUU|) look simpler 

+ -^)f + + M)g = 0, 

dp sm p cos p 

i± - - (^ - M) / = . (102) 

dp sm p cos p 

Summing and subtracting two last equations, we get 

^if + 9) + ^{f-g)-^{f-9)+Mif + g)=0, 
dp sm p cos p 

^if-g) + ^if + 9) + ^if + 9)-Mif-g)=0. (103) 
dp sm p cos p 

Introducing two new functions 

f + g = e-^P/\F + G), f-g = e+'^/^F - G) , (104) 
one translates (|103|) into 

^e-'P'^{F + G) + -^e+'P'^iF - G) —e+^P''^{F - G) + Me-'P/^(F + G) = , 

dp sin p cos p 

—e+'P/^(F - G) + -^e-'P/^(F + G) + -^e-'P/'^(F + G) - Me+'P/'^(F - G) = , 
dp sin /9 cos p 



or 



(F + G) - ^ (F + G) + (cos /) + i sin />) (F - G) 
dp 2 sm p 

^ -(cosp + isinp)(F-G)+M(F + G) = 0, 



cos p 



^{F - G) + ^(F - G) + -^(cosp - i sin p)(F + G) 
dp 2 sm p 

+— ^(cosp-isinp)(F + G) -M(F- G) = . 
cos p 

Now summing and subtracting two last equations, we get 

, d cos p sin p , , i . „ 

+ — ^ - ie ^ F+ e + M-iz.--G = 0, 

dp sm p cos p 2 

(A_,^ + ,,!i!i^)G + (-e + M + ..-l)F = 0. (105) 
dp sin p cos p 2 



This system can be simplified by substitutions 



F(p) = sin"'' p cos"^' p v?(p) , G(p) = sin+'^ p cos+^' p r(p) , 

d i 
sin"^'' pcos"^*' p— (p + (e + M - ii^ - -)r = , 
dp 2 

d i 
sin+2'' pcos+2" p— r +(-e + M + iu--)(f = . (106) 
dp 2 
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Let us specify a 2-order differential equations for ip 

<^ , ^ cosp ^. sinp, d . .n . ,„N2n 

+ {-2v + 2%e -)— + (e - ivf - (M - i lf\ = . 

sinp cos p dp 

Turning back to the variable r = sin p 

dr"^ r dr 

and translating the equation to the variable z = r'^ one obtains 

^(1 -z)^ + [{^-u)-{l-u- ie)z]^ - i[(M - z/2f - (6 - = . (107) 

The later is of hypergeometric type with parameters defined by 

7=^-i/, a + l3 = -u-ie, aP = ^[{M - i/2f - (e - iuf] ; 

from whence it follows 

_ -u-ie + {iM + 1/2) -u - ie - {iM + 1/2) 

2 ' ^" 2 ■ 

Thus, solution are are given by (remembering, v = j + 1/2) 

-j - i(e + M) -i - 1 - i(e - M) 

2 ' ^ = 2 ' 

7 = -i = -l/2,-3/2,..., (^(r) = F(a,^,7,-2 = 0. (108) 

The main conclusions about penetration of the spin 1/2 particle through de Sitter horizon 
remain the same: they coincides with results for S = 0, 1 fields. 
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